OPEN Q ACCESS Freely available online 



&PLOS 



ONE 



Kinetic Titration Series with Biolayer Interferometry 



Daniel Frenzel 1 , Dieter Willbold 1 



1 Forschungszentrum Julich, ICS-6 Structural Biochemistry, Julich, Germany, 2 Heinrich-Heine-Universitat Dusseldorf, Institut fur Physikalische Biologie, Dusseldorf, 
Germany 



Abstract 

Biolayer interferometry is a method to analyze protein interactions in real-time. In this study, we illustrate the usefulness to 
quantitatively analyze high affinity protein ligand interactions employing a kinetic titration series for characterizing the 
interactions between two pairs of interaction patterns, in particular immunoglobulin G and protein G B1 as well as scFv IC1 6 
and amyloid beta (1-42). Kinetic titration series are commonly used in surface plasmon resonance and involve sequential 
injections of analyte over a desired concentration range on a single ligand coated sensor chip without waiting for complete 
dissociation between the injections. We show that applying this method to biolayer interferometry is straightforward and i) 
circumvents problems in data evaluation caused by unavoidable sensor differences, ii) saves resources and iii) increases 
throughput if screening a multitude of different analyte/ligand combinations. 
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Introduction 

Surface plasmon resonance (SPR) is widely used to study 
analyte/ligand interactions in real-time [1]. For SPR analysis, one 
interactant ("ligand") is immobilized onto the sensor surface while 
the other interactant ("analyte") is passed over this surface by 
continuous flow. The standard assay requires complete removal of 
the analyte ("regeneration") between each measurement cycle to 
avoid any residual analyte molecules blocking potential binding 
sites on the surface. This procedure is termed "multi-cycle 
kinetics" [2], because it consists of several cycles of alternating 
analyte injections and surface regeneration steps. The regenera- 
tion process often requires conditions that can inactivate the 
immobilized ligand irreversibly [3]. To avoid such potentially 
detrimental surface regeneration steps, a technique has been 
developed that allows several concentrations to be applied within a 
single cycle without the requirement of complete surface 
regeneration steps following each injection. These so called 
"kinetic titration series" or "single cycle kinetics" were found to 
be virtually as precise as classical "multi cycle kinetics" [4] and can 
be used as an alternative option if regeneration is not practical 
[5,6] . Later kinetic titration series were adopted by GE with minor 
modifications. and renamed as "single-cycle kinetics" (Biacore 
manual). 

More recently, biolayer interferometry (BLI) has become an 
alternative method to SPR. One advantage of BLI is that the 
number of sensors can be scaled up easily without making the 
system more error-prone or complex. Theoretically, there is no 
need to regenerate single sensors, because duplicates of the surface 
can be easily created by immobilizing an equal amount of ligand 
on additional sensors. However this approach has some disadvan- 
tages. It is not possible to achieve identical ligand coatings of 
multiple sensors. Such an approach also increases sensor 



consumption and it is not guaranteed that each sensor has an 
equal performance in later measurements. 

We strove to overcome these disadvantages of BLI as compared 
with SPR by exploring whether application of a kinetic titration 
series in BLI experiments is feasible and accurate. 

Materials and Methods 

Preparation of protein G B1 (GB1) and scFv IC16 

Purification of GB 1 was done by standard Immunoglobulin G 
(IgG) affinity purification after expression in E. call with pGEV2- 
GB1 [7] (see SI: "Preparation of protein G Bl") for details. 
Purification of scFv IC16 was done as described in Frenzel et al. 
[5]. 

Immobilization of ligands via amine coupling on (AR2G) 
biosensors 

The 40 LlM GB1 solution in PBS was diluted in 10 niM sodium 
acetate buffer pH 4.0 to a final concentration of 20 U.M (binding 
buffer). The sample sensors were pre-incubated in ddH a O for 
10 min, activated in a 1:1 mixture of 0.1 M N-Hydroxysuccini- 
mide (NHS)/0.4 M l-Ethyl-3-(3-dimethylaminopropyl)-carbodii- 
mide (EDC) for 800 s and incubated in binding buffer for 900 s. 
The reference sensors were activated in the same way, but not 
treated with binding buffer. All sensors were blocked with 1 M 
ethanolamine for 180 s and stored in ddF^O before further usage. 
All steps were performed at 20°C with an agitation speed of 
1000 rpm. For method definition and execution, the Data 
Acquisition software 7.1.0.92 from ForteBio was used. 
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Immobilization of C-terminally biotinylated AG(1-42) via 
amine coupling on Super SA-biosensors 

C-terminally biotinylated AF3(l-42) (EUROGENTEC) was 
dissolved in 100% HFIP and incubated at RT overnight. The 
stock solution was divided in 26.5 fXg aliquots. HFIP was removed 
by evaporation in a Concentrator 5301 (Eppendorf). A(5 was 
freshly solubilized in 550 sodium phosphate buffer pH 7.4 
(10 mM; yielding a 10 |lM A|3 solution). To separate the 
monomers from bigger particles, they were subjected to a density 
gradient centrifugation (DGC) as described in Frenzel et al. [5]. 
After centrifugation, fraction one (140 JLll) was used for immobi- 
lization of AB(l-42) monomers via standard strep tavidin-biotin- 
coupling procedure with Super SA-sensors (ForteBio). It was 
planned to immobilize 0.15 nm, 0.25 nm and 0.75 nm of ligand 
on eight sensors respectively. Further eight sensors were used as 
reference and remained in phosphate buffer. The finally achieved 
layer thickness of all sensors is summarized in Tab. 1. 

Parallel sensor kinetics of GB1 with biolayer 
interferometry 

Kinetic titration series were performed in the interaction buffer 
(PBS with 0.05% Polysorbat 20). 5 mg/ml IgG (ID: ABIN376828; 
Antibodies-Online) was diluted in interaction buffer to 0.5 (J.M 
and further diluted four times with a dilution factor of two. To 
measure the interaction between IgG and GB1, the association 
and dissociation times were 360 and 600 s, respectively, for every 
analyte concentration. In total, five sensors were used to measure 
five different analyte concentrations in parallel, while one sensor 
was used to measure the buffer reference. Additional six sensors 
were used as sensor reference. All steps were performed at 25°C 
with an agitation speed of 1 000 rpm. Sensorgrams were measured 
on an Octet Red96 (ForteBio) and double referenced against the 
buffer reference signal and the reference sensor signals using the 
Data Analysis software 7.1.0.36 (ForteBio). The double referenced 
sensorgrams were exported into the BiaEvaluation 4. 1 compatible 
"csv"-format by a python script (see SI: "Scripts"). The 
sensorgrams obtained with the concentrations: 0.5, 0.25, 0.125, 
0.0625 and 0.03125 [LM were fitted with the BiaEvaluation 
software 4. 1 from Biacore using a 1 : 1 binding model that included 
an Rl-term. 

Parallel sensor kinetics of scFv IC16 with biolayer 
interferometry 

Kinetic titration series were performed in the interaction buffer 
(PBS with 0.5% Polysorbat 20, 0.1% BSA and 10% NSB reducer 
from GE Healthcare). 2.4 |0,M scFv IC16 was diluted four times 
with a dilution factor of two. To measure the interaction between 
AB(l-42) and scFv IC16, the association and dissociation times 
were 270 and 90 s, respectively, for every analyte concentration. 
Further steps are comparable with "Parallel sensor kinetics of GB1 
with biolayer interferometry" (Data Analysis software: 8.0.0.35). 
The sensorgrams with the concentrations: 0.24, 0.12, 0.06, 0.03 
and 0.015 |J,M were fitted with the BiaEvaluation software 4.1 
from Biacore using a 1 : 1 binding model without Rl-term. 

Kinetic titration series of GB1 with biolayer 
interferometry 

To measure the affinity between IgG and GB 1 , the association and 
dissociation phases were recorded for 360 and 240 s, respectively, for 
every analyte concentration (same concentrations as described in: 
"Parallel sensor kinetics of GB 1 with biolayer interferometry"). Four 
sensors recorded the kinetic titration series, whereas one sensor 
recorded the buffer reference signal. Additional five sensors were used 
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Figure 1. Evaluation of (A) kinetic titration series and (B) parallel sensor kinetics with rabbit IgG binding to GB1 in BLI. The 

sensorgrams show the interaction of IgG (analyte) with GB1 (ligand). Applied analyte concentrations were: 0.5, 0.25, 0.125, 0.0625 and 0.03125 |iM. 
The fits are indicated by the red lines, whereas the sensorgrams are shown in black (A) and blue (B). The residuals of the fits are plotted below the 
respective sensorgram. All other experiments are shown in File SI. 
doi:1 0.1 371 /journal.pone.01 06882.g001 



as sensor reference. All steps were performed at 25°C with an 
agitation speed of 1000 rpm. The sensorgrams were double 
referenced against the buffer reference signal and the empty sensors 
by the Data Analysis software 7.1.0.36 (ForteBio). The double 
referenced signals of each association and dissociation phase were 
combined and exported into a BiaEvaluation 4. 1 compatible "csv"- 
format using a python script (SI: "Scripts"). The sensorgrams were 
fitted with the BiaEvaluation software 4. 1 from Biacore with a 1 : 1 
kinetic titration series model that included an Rl-term [4] . 



Kinetic titration series of scFv IC16 with biolayer 
interferometry 

To measure the affinity between IgG and GB 1 , the association 
and dissociation phases were recorded for 270 and 90 s, 
respectively, for every analyte concentration (same concentrations 
as described in "Parallel sensor kinetics of scFv IC16 with biolayer 
interferometry"). Five sensors recorded the kinetic titration series, 
whereas one sensor recorded the buffer reference signal and six 
sensors were used as sensor reference. The other steps are 
comparable with section "Kinetic titration series of GB1 with 
biolayer interferometry" (Data Analysis software: 8.0.0.35). The 
sensorgrams with the concentrations: 0.24, 0.12, 0.06, 0.03 and 
0.015 |XM were fitted with the BiaEvaluation software 4.1 from 



PLOS ONE | www.plosone.org 



3 



September 2014 | Volume 9 | Issue 9 | e1 06882 



Kinetic Titration Series with BLI 



A 




200 400 600 800 1000 1200 1400 1600 1800 
time in s 



150 200 
time in s 




S -0.05 



'V, 



B 




time in 5 





0.10 


E 


c 




c 


0.05 


uals 


0.00 






s 


-0.05 



800 1000 1200 1400,' 1600 1800 





0.10 


E 




c 




c 


0.05 


uals 


0.00 






S 


-0.05 



IS*. 



150' ' 200 250 "§W W ^50" 




800 1000 
time in s 



150 200 
time in s 



0.20 

= 0.10 

E 

la 0.00 
n i 

§ -0.10 0 

£> -0.20 



800 1000 1200 





0.20 


E 




c 


0.10 


c 






0.00 


~nj 




T! 


-0.10 


o 


-0.20 




250 ~3Kr*"~*"350 



Figure 2. Comaprison of kinetic titration series (A-C) and parallel sensor kinetics (D-F) with scFv IC16 binding to A|i(1-42) in BLI. 

The sensorgrams show the interaction of scFv IC16 (analyte) with C-terminally biotinylated Afi(1-42) (ligand). The amount of ligand was increased 
from 0.1 3 nm (A, D), 0.41 nm (B, E) and 1 .01 nm (C, F). Applied analyte concentrations were: 2.4, 1 .2, 0.6, 0.3 and 0.1 5 uM. The fits are indicated by the 
red lines, whereas the sensorgrams are shown in blue. Each kinetic titration series was reproduced five times. The residuals of the fits are plotted 
below the respective sensorgram. 
doi:1 0.1 371 /journal.pone.01 06882.g002 
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Biacore with a 1:1 kinetic titration series model without Rl-term 
[4]. 

Results 

Ligand immobilization 

In order to compare the practicality and efficiency between 
multi cycle kinetics and kinetic titration series using BLI, we used 
the well-studied interaction of GB1 with IgG. The reported 
dissociation coefficients (K D ) of GB 1 to the constant (Fc) region of 
IgGs are -0.1 |lM (human IgG) and -0.77 nJVI (rabbit IgG) 
[8,9]. The interaction of the scFv IC16 with Ap(l— 42) was used as 
an additional example system. ScFvs show in comparison to IgGs 
no avidity. The scFv IC16 is directed against the N-terminus of 
AP(l-42) and with SPR, a reliable K D value of 0.76 uM for C- 
terminally biotinylated A(5(l-42) monomers was already estimated 
[5]. 

GB1 was immobilized onto AR2G sensors via amine coupling. 
We found that at the end of the multiple immobilization 
procedures the amount of immobilized protein differed for each 
sensor (Tab. 1). The mean and standard deviation was 1.24 nm 
and 0.20 nm. Thus, homogeneous immobilization of protein to 
the surface was not possible, because the on-rates of the sensors 
seem to deviate from each other. We estimated fewer deviations 
with Streptavidin-Biotin coupling on "Super SA" sensor tips, 
especially at higher layer thicknesses. The sensors with 1.01 nm 
ligand had a standard deviation of 0.05 nm, which corresponds to 
4.9%, based on the mean (Tab. 1A). However, with declining 
amount of ligand, the ratio of standard deviation to immobilized 
ligand grew (see Tab. IB: 4.8% for 0.41 nm, SD: 0.02 nm and 
Tab. 1C: 15.4% for 0.13 nm, SD: 0.02 nm). 

Binding kinetics with parallel sensor kinetics 

Parallel sensor kinetics is thought to be more precise than a 
kinetic titration series, because no secondary processes (like 
.£ dissociation of previously bound analyte from the surface) impair 

% the measurements. However, other effects like inhomogeneous 

c coating and differences in sensitivity of single sensors are expected 

^ to compromise these precision advantages. 

To obtain data from parallel sensor kinetics, five sensor pairs 
o were applied to record the signals received from applying five 

£ different analyte concentrations. One sensor pair was applied to 

£ record the buffer reference signal. Measurements with IgG and 

£ GB1 were fitted globally with Rl-term (see: Fig. 1), whereas 

i= measurements with scFv IC16 and AfS(l-42) were fitted without 

~ Rl-term (see Fig. 2). The obtained K D s for GB1 were 0.16 JiM 

2 and 0.25 (XM. The % 2 values were 5.04xl0~ 6 nm 2 and 

0 2.38 xlO - nm . The term % 2 gives a measure for the accuracy 
of the fitting [10]. It represents the averaged, squared residual per 

5 data point. In our case, % 2 was below the squared sensor noise 

•- (—0.008 nm), which is a quality indicator of a fit. With the scFv- 

S system, we obtained 0.18 \xM/% 2 : 7.17xl0~ 5 for 0.15 nm ligand, 

% g 0.59 uM/x 2 : 2.25 xlO" 4 for 0.41 nm ligand and 0.43 uM/x 2 : 

JjS 9.69 xlO" 4 for 1.01 nm ligand (experiment was reproduced: data 

g gg on request). Remarkably the K D s spread and just the best fit (see 

TO O 

1 °. 

"5 o 



Fig. 2: E) is close to the expected affinity range [5]. 



Binding kinetics with kinetic titration series 

ro £ Sensorerams from kinetic titration series were recorded for each 

£ amount of ligand from five (scFv IC16 and AB1-42) and seven 

^ ^ (IgG and GB1) sensor pairs (one sensor with ligand and one 

without) that were each subsequently applied to five different 
analyte concentrations within one titration series. Additionally one 
sensor pair was applied to record the buffer reference signals. Like 
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for parallel sensor kinetics, measurements with GB1 were fitted 
with Rl-term, whereas measurements with scFv IC16 were fitted 
without Rl-term. Each sensor pair was fitted separately and each 
fit was used to calculate the mean values for the rate constants A a 
and k d> as well as the dissociation constant K D . The estimated 
mean K D of the interaction of IgG and GB1 was 0.16 |iM (7 
replicates, standard deviation: 0.08 |0,M) and the obtained mean % 2 
was 4.92xl0~ 6 nm 2 (standard deviation: 2. 15 x 10" b nm 2 ). By 
taking into account the standard deviation, the % 2 is virtually 
identical to the % 2 obtained for parallel sensor kinetics. Compar- 
ison of the on-rates (& a ), off-rates (k^j and binding constants (^d) 
for IgG and GB1 showed that both methods give near identical 
values (Tab. 2). Both, the on-rates and the off-rates lie within the 
same range of 10 4 1/Ms and 10 3 1/s respectively (Fig. 1). 

With the scFv-system we estimated the K D s: 0.54 mJVI/% 2 : 
3.87xl0~ 5 for 0.15 nm ligand, 0.64 uM/% 2 : 6.80xl0~ 5 for 
0.41 nm ligand and 0.71 [lM/% 2 : 2.44xl0" 4 for 1.01 nm ligand 
(experiment was reproduced with three sensor pairs: data on 
request). It is obvious, that for every amount of ligand, the x 2 -term 
is at least several times smaller in comparison to parallel sensor 
kinetics (see Fig. 2: A/B/C vs Fig. 2: D/E/F) and all the estimated 
K D s (see Tab. 2 and Tab. SI in File SI) are very close to the value 
estimated by SPR [5]. 

A known problem of using the % 2 rating is that this method is 
based on averaging. Local deviations can hardly be evaluated with 
% 2 alone. If one fit with a higher / 2 and one with a smaller % 2 are 
compared with each other, it is possible that the fit with the smaller 
% 2 has higher local deviations from the sensorgram. In contrast, 
the fit with the higher x 2 could proceed completely within the 
noise pattern. In this example, the fit with the higher % 2 could be 
the more accurate description of the sensorgram. In our case, the 
fit curves of the kinetic titration series are frequently within the 
range of the sensor noise, whereas the fit of the parallel sensor 
kinetics is outside the noise range at certain time points (Fig. 1 A/ 
Fig. SI in File SI: t 360 „, t 960 s , t 1560 s , t 2 i 60 s , Fig. IB/Fig. S2 in File 
SI: to_ 6 os and t 36 o-4oo s and t 2 7eo s ; Fig. 2). This illustrates that the 
kinetic titration series can yield more reliable fits, because the 
affinity differences of single sensors are omitted since only one 
sensor pair is used per interaction study. This is most obvious if the 
sensorgrams are fitted without Rl-term or local R max as a linear 
correction mechanism (Fig. 2: A-F and Fig. S3 in File SI). 

Discussion 

For our example system rabbit IgG and GB1, kinetic titration 
series and parallel sensor kinetics provided near identical results 
with regard to the on/ off-rates as well as the K D values and are in 
accord to published data. The mean % 2 values (considering all 
sensor data of either the kinetic titration series or parallel sensor 
kinetics) were nearly identical in both methods. The other example 
system, utilizing the interaction of scFv IC16 and A|3(l-42) was 
fitted without Rl-term. The result was showing a clear advantage 
of the kinetic titration series in respect to x 2 -values the reliability of 
the estimated K D at every immobilization level. However, it is not 
advisable to rate fits based only on the % 2 . This is why we consider 
local deviations of the fits from the sensorgrams as another marker 
for the quality of the fit. With regard to this point, fits of parallel 
sensor kinetics have stronger local deviations from the sensorgrams 
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